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ABSTRACT 

We present Subaru/HiCIAO //-band high-contrast images of the debris disk around HIP 79977, whose pres- 
ence was recently inferred from an infrared excess. Our images resolve the disk for the first time, allowing 
characterization of its shape, size, and dust grain properties. We use angular differential imaging (ADI) to 
reveal the disk geometry in unpolarized light out to a radius of ~2", as well as polarized differential imag- 
ing (PDI) to measure the degree of scattering polarization out to ^1"5. In order to strike a favorable bal- 
ance between suppression of the stellar halo and conservation of disk flux, we explore the application of 
principal component analysis (PCA) to both ADI and reference star subtraction. This allows accurate for- 
ward modeling of the effects of data reduction on simulated disk images, and thus direct comparison with 
the imaged disk. The resulting best-fit values and well-fitting intervals for the model parameters are a sur- 
face brightness power-law slope of S out — —3.2 [—3.6,-2.9], an inclination of i = 84° [81°, 86°], a high 
Henyey-Greenstein forward-scattering parameter of g = 0.45 [0.35,0.60], and a non-significant disk-star off- 
set of u = 3.0 [—1.5,7.5] AU = 24 [—13,61] mas along the line of nodes. Furthermore, the tangential linear 
polarization along the disk rises from ^10% at 0"5 to ~45% at 1"5. These measurements paint a consistent 
picture of a disk of dust grains produced by collisional cascades and blown out to larger radii by stellar radiation 
pressure. 

Subject headings: circumstellar matter — planetary systems — techniques: high angular resolution — stars: 
individual (HIP 79977) 



1. INTRODUCTION 

Most debris disk systems are second-generation disks pro- 
duced by the collisional destruction of planetesimals, al- 
though some of the youngest systems (^10 Myr) might 
also include primordial star formation material (e.g., Wyatt 
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2008). While space-based inf rared (IR) surveys ha ve iden- 
tified many debris disks (e.g., Auman net al.l 1 1984b . though 
few have been spatially resolved at optical or IR wave- 
lengths. High-resolution images have revealed morpholog- 
ical disk sub-struct ures such as warps, geometrical offsets , 
and clumps (e.g iKalas & Jewittl fl995b iHeap et alJ l200(il : 
lKalasetalJ[2005l [Lagrange et al. 2009), possibly caused by 
gravitational perturbations from exoplanets. In particular, 
identifying new debris disk systems with an edge-on incli- 
nation is critical to efforts to study the origin and compo- 
sition of their small, second-generation gaseous components 
(Roberge et al.ll2006l) . 

HIP 79977 (HP 146897) is a 5-10 Myr old (Ide Zeeuw et alJ 
[19991 iPecaut eTaTJ 12011 iSong et alJ I20T2I) F 2/3V member 
of Up per Scorpius at a distance of 123^4 pc (Ivan Leeuwenl 
120071) which has been identified as a debris disk system with 
a fractional IR luminos ity Lir/L* = 5.9 x 10~ 3 dChen et all 
l2006t IChen et aTJ 120 1 ll) . A simple single-temperature fit to 
its IR excess led IChen et al.l d201 lb to suggest a source ring 
radius of at least 40 AU, with a characteristic temperature of 
89 K. Here we present high-contrast imaging of HIP 79977 
which resolves this debris disk in scattered light for the first 
time. 

2. OBSERVATIONS 

As part of the SEEDS survey (Strategic Exploration o f Ex- 
oplanets and Disks with Subaru/HiCIAO: lTamurall2009l) . we 
obtained two epochs of high- contrast imaging on HIP 79977 
with the HiCIAO instrument dHodapp et al.| [2008) on Subaru 
Telescope. 

The first dataset was taken in //-band (1.6 /xm) on 2012 
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FIG. 1. — Subaru HiCIAO //-band images of the HIP 79977 debris disk, (a-c) Logarithmic intensity images based on the May 2012 data as produced by the 
reduction methods described in the body text, (a) Classical ADI. (b) Conservative LOCI (Ng = 0.75, Na = 10,000). (c) PCA-assisted ADI (5 eigenmodes). (d) 
Signal-to-noise (S/N) map for image c in linear stretch. The area outlined in yellow was excluded for noise calculation. SI marks a likely background star, (e) 
Cleaned polarized intensity image from the luly 2012 data, in logarithmic stretch, (f) Linear polarization orientation angle. The mean orientation of the disk flux 
averaged in 15-pixel boxes along the midplane is 25.7° ± 3.5°, consistent with the expected tangential orientation of 24.0°. All images were convolved with a 
5-pixel (« 1 FWHM) diameter circular aperture. 



May 12, comprising 119 frames of 20 seconds for a total 
exposure time of 39.7 min, and spanning 19.6° of field rota- 
tion. The image rotator was operated in pupil-tra cking mode 
to en able angular differential imaging (ADI; Mar ois et al.l 
2006). The plate scale was 9.50 ±0.02 mas per pixel, and 
the field of view 20"x20 ". The AO 188 adaptive optics sys- 
tem dMinowa et al.120 10) provided a stable point-spread func- 
tion (PSF) with a FWHM of 6.7 pixels = 64 mas under good 
weather conditions (~0"6 DIMM seeing). HIP 79977 sat- 
urated out to a radius of 22 pixels = 0"2; no coronagraph 
was used. After co rrection for flatfield and field distortion 
dSuzuki et al.l 120 10). the images were registered by fitting a 
Moffat profile to the PSF halo, for an estimated registration 
accuracy of 0.5 pixels = 5 mas. 

The second dataset was taken in two-channel polarized dif- 
ferential imaging (PDI) mode on 2012 July 7 under the same 
filter, plate scale, and saturation conditions, with excellent 
seeing (0"4). It comprises 64 frames of 30 seconds for a to- 
tal exposure time of 32 min, and spans 17.3° of field rotation. 
The field of view was 10"x20". The exposures are organized 
into serial batches of four frames in which HiCIAO's half- 
wave plate cycled through the position angles (0°, 45°, 22.5°, 
67.5°), measuring the Stokes linear polarization parameters 
(7+ Q,I— Q,I + U,I— U), respectively. This allows for cor- 
rection of non-common path aberrations by double-difference 
polarimetry on short time scales. Each frame was registered 
by Moffat fitting. While the intensity component of the PDI 
data set can also be reduced with ADI, the resulting quality 
is inferior to that of the dedicated May 2012 ADI data. We 
therefore use the May 2012 data to constrain the full-intensity 
appearance of the HIP 79977 disk, and the July 2012 data to 
obtain the matching polarized-intensity image. 



3. DATA REDUCTION 
3.1. Angular differential imaging 

While the API technique with the LOCI algorithm 
dLafreniere et al. 2007|) is commonly used to improve high- 
contrast sensitivity for detecting extra solar planets (e.g., 
Maroi set al.l 120101: [L agrange et al.l 120101) and brown dwarf 
companions (e.g.. IThalmann et aD l2009). we recently demon- 
strated its use for detecting faint circumst ellar disks hid- 
den in the speckle halo of their host stars (iThalmann et al.l 
IMotlBuenzhetal.ll2010tlThalmannetalJl2011l) . This "con- 
servative LOC I " tech n ique has sin c e been widely adopted 
dRodigas et all 120121; ICurrie et al.l 120121: iBoccaletti etaTI 
120 12t Lagrange et al.ll2012l) . 

The benefit of powerful speckle suppression with ADI 
comes at the price of partial lo ss of disk flux and morpholog- 
ical integrity dMilli et al.l 12012b . However, careful forward- 
modeling of these effects can recover the physical disk mor- 
phology, as we demon strated on the HR 4796 A debris disk 
(Thal mann et al.ll201 ll) . Clarification is needed regarding re- 
cent discussions that attributed some results based on this 
techniq ue to ADI artifacts rather than to ast rophysical pro- 
cesses (lLagrange et al.ll20Tl iMilli et al.ll20l2l) . Since ADI 
works in concentric annuli, it cannot generate spurious sig- 
nals at a radius where no disk flux is present. The pres- 
ence of "streamer"-like flux distributions therefore does imply 
the presence of physical emission at those radii, even though 
the streamer morphology itself is known to be carved from a 
p hysically smooth disk b y the ADI process. The conclusions 
in IThalmann et al.l (1201 ll) therefore remain valid, and we con- 
tinue to use this methodology in this work. 

For HIP 79977, we started with classical ADI dMarois et all 
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2006, no frame selection) and conservative LOCI (Ng = 0.75, 
Na = 10,000). Furthermore, we adopted the application of 
principal component analysis (PCA) to API as a robust an d 
deterministic alternative to L OCI (|Soummer et al.l 120121) . 
Closely following the recipe of Amara & Ouanz (2012|), we 
performed PCA on the stack of centered, pupil-stabilized 
frames after subtracting the mean image and masking 
the saturated central region (r < 22 pixels = 0"2). We 
then projected each frame onto the first n modes of the 
orthonormal set of Karhunen-Loeve eigenimages delivered 
by PCA, and subtracted those projections from the frame. 
The stack of frames was then derotated and collapsed to 
yield the final image. We hereafter refer to this method as 
PCA-ADI. Through visual inspection of the resulting images, 
we found that n = 5 provides an excellent trade-off between 
suppression of the stellar PSF halo and conservation of disk 
flux. Note that PCA-ADI with n = is identical to classical 
ADI. To ensure linearity of our data reduction process, we 
use mean-based rath er than median-bas ed frame combination 
in all techniques (cf. Brandt et al. 2 0121) . 

Finally, we explored the use of PCA for PSF subtraction 
with reference stars. However, we found that for HiCIAO H- 
band data, the speckle halo varied too much between the sci- 
ence and the reference targets to reach competitive contrasts. 

Scattered light from the debris disk around HIP 79977 is 
clearly detected in the May 2012 data set with all ADI-based 
techniques, as shown in Figure QJ-c. The ripple pattern vis- 
ible in the first two images is a spurious residual of Sub- 
aru's spider diffraction pattern. Only the PCA-ADI technique 
proves effective at removing this artifact and revealing the 
smooth slope of the debris disk. We therefore choose it as 
the benchmark technique for further analysis, as discussed in 
Section |4] Figure [TJl shows the signal-to-noise (S/N) map of 
the PCA-ADI image, calculated as the standard deviation in 
concentric annuli after smoothing with a 5-pixel diameter cir- 
cular aperture. A box containing the disk flux was masked for 
noise calculation. 

The disk presents itself almost edge-on, and thus appears 
as a ribbon in the ADI images. Most of the flux is concen- 
trated to the South of the projected disks' major axis, identi- 
fying this as the disk's likely near side and indicating strong 
anisotropic forward-scattering. The disk flux remains statisti- 
cally significant from ~0"3 (« 40 AU) out to -2" (« 250 AU 
in projection). 

3.2. Polarized differential imaging 

We reduced the July 2012 polarized differential imaging 
data using the sta ndard procedure for SEEDS polarimetry (cf . 
iTanii et all 120 121) . It includes double-difference polarimetry 
using four half-wave plate position angles, as well as mod- 
eling and corre ction of instrumental polarization based on 
Uoos et all (12008b . The intensity PSF was fit to the Q and 
U images and subtracted to remove contributions from a po- 
larized stellar PSF. The resulting polarized intensity image 
P = V (Q 2 + U 2 ) was found to contain spurious structures left 
behind by the subtraction of imperfectly aligned PSFs, which 
we removed by projecting spatial derivatives of the intensity 
PSF to the Q and U images, and subtracting them. Finally, we 
noticed that RMS of the positive and negative residual speckle 
noise in the Q and U images added up to a positive radial halo 
in the P image. We measured this noise profile at position an- 
gles that do not show disk flux, and subtracted it in quadrature 
from the entire image. 




FIG. 2.— Modeling of the HIP 79977 debris disk, (a) Simulated logarith- 
mic scattered-light image of the best-fit model disk, (b) The same after ap- 
plying PCA-ADI data reduction (linear stretch to show oversubtraction). (c) 
Noise-normalized residual map of the PCA-ADI reduction of the May 2012 
data, at a linear stretch of [— 2.5o~,5o~]. The evaluation region is marked 
with a yellow outline, (d) The same calculated after subtraction of the model 
disk treated with PCA-ADI. The disk flux is effectively removed from the 
data (residual Xn\m = 1095 = l.lGAfcata). An unconfirmed point-like signal at 
4.6 o~ significance is highlighted with an arrow. For visual clarity, no binning 
has been applied for panels (c,d). 

The resulting P image shows the debris disk out to a sep- 
aration of r « 1"5, where the polarized flux level vanishes 
into the background (FigureQJ). The polarization vectors are 
tangentially oriented, consistent with scattering polarization 
(Figure QJ). No brightness contrast between the near and far 
side of the disk is evident, in accordance with the expectation 
that the excess scattered light from forward scattering on large 
grains is unpolarized. 

4. ANALYSIS 

4.1. Disk modeling 

In order to extract physical disk information from our PCA- 
ADI image, we generate simulated scattered-light images of 
model disks using the GRaTer code (Augereau et al. 1999). 
The main free parameters of the model are the semimajor axis 
of the source ring ao, the inner and outer power-law slopes 
OCimOCout of the radial density distribution, the disk-star off- 
set //. the inclination angle i, the position angle <p, and the 
Henyey-Greenstein parameter g characterizing the anisotropy 
of the scatt ering behavior of the d ust grains. For definitions 
we refer to Au gereau et al.l (11999b . Note that the measurable 
outer slope of the disks radial surface brightness profile is not 
Ooub but Som = OJcmt + fi — 2, where j3 represents the radial 
shape of the scale height distribution. While this relationship 
breaks down for purely edge-on disks, it holds for the radial 
profile for our models with i < 88°. We fix j3 = 1 and hence- 
forth use S ou t t0 characterize the disks outer slope. Since we 
do not resolve the gap inside the source ring, we adopt the 
semimajor axis ao = 40 AU (~ 0"3) derived from infrared ex- 
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TABLE 1 

Summary of Properties and Results for HIP 79977 
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FIG. 3. — Constraints on the disk model parameters S out (outer slope), i 
(inclination), and g (Henyey-Greenstein parameter). The contours delimit 
the parameter space of well-fitting solutions (% 2 • 
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cess bv lChen et al] rf20ll . 

As a first step, we determine the position angle (j) of the 
projected disk's axis of symmetry. We rotate the PCA-ADI 
image from its North-up, East-left orientation by a given an- 
gle, mirror this image about the y axis, and subtract it from an 
unmirrored copy, adjusting the angle to minimize the subtrac- 
tion residuals. To quantify the residuals, we first bin the resid- 
ual image by a factor of 5 (« 1 FWHM) in both dimensions, 
and define an evaluation region of A'data = 1025 binned pixels 
framing the disk flux (Fig. |2j. We derive a noise profile as 
the standard deviation of pixel values in concentric annuli in 
the residual image, excluding the evaluation region. We then 
divide the residual image by this noise profile and calculate 
X 2 as the sum of squares of all pixel values in the evaluation 
region. 

Our data reduction does not yield a residual map in which 
the pixels have fully independent normal errors, even after 
binning. Therefore, we do not expect the usual % 2 thresh- 
olds to produce reliable confidence intervals, and use a much 
more conservative % 2 threshold of \/2A^jata- Rather than de- 
termining confidence intervals about a best-fit model, we seek 
to determine the family of models consistent with our images. 

For the position angles this yields best-fit value of = 24.0° 
(counter-clockwise from North) and a well-fitting range of 
[23 .7°, 24.3°]. This includes HiCIAO's True North offse t of 
0.35° ±0.02° calibrated using lvan der Marel et all (120071) . 

Next, we attempt to match the observed disk morphology 
by calculating a grid of model disks with i = [77°... 88°], 
Sout = h2.8 ... - 3.7], and g = [0.00. . .0.80]. Rather than 
implanting the mode ls into the raw data as we did in 
Thalm ann et alj d201 ll) . we exploit the deterministic nature of 
PCA-ADI to forward-model the ex act effect of our data re duc- 
tion process on the model disk (cf. lSoummer et all2012l) . We 
subtract a scaled version of this processed disk image from 
the data, choosing the scale factor to minimize residual % 2 as 
defined above. This yields best- fit model parameters and well- 
fitting parameter ranges of = 84° [81°, 86°], 5 0U t = -3.2 
[-3.6,-2.9], g = 0.45 [0.35,0.60], and a ratio of total disk 



Fixed disk model parameters 


adopted 
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Source ring semimajor axis ciq (AU, ")''' 


40 ~ 0"33 


< 60 ~ 0"49 


Inner brightness distribution slope a m 


20 




Scale height radial shape parameter j8 


1 




Optimized disk model parameters 


best fit 


well-fitting 


Minor axis position angle (°) 


24.0 


[23.7,24.3] 


Inclination i (°) 


84 


[81,86] 


Henyey-Greenstein parameter g 


0.45 


[0.35,0.60] 


Disk-star offset u (AU)< 2) 


3.0 


[-1.5, +7.5] 


Eccentricity e = abs(u/oo) 


0.06 


[0,0.16] 


Outer brightness distribution slope S out = 


Oou, + ft - 2: 




— model disk, true slope 


-3.2 


[-3.6,-2.9] 


— model disk, after PCA-ADI 


-2.6 


[-2.9,-2.3] 


— measured, after PCA-ADI 


-2.7 


[-2.8,-2.6] 


Disk/star flux contrast, //-band 


0.0033 


[0.0029,0.0038] 


Degree of scattering polarization 


measured 


1 o" interval 


— at0"5~62AU(%) 


10 


[5,20] 


— at l!'0ss 123AU(%) 


35 


[20,45] 


— at l'.'5~ 185AU(%) 


45 


[30,60] 



NOTES. The given ranges for parameters ;', g, S out , e, and u represent the min- 
imum and maximum values of the well-fitting family of disk models defined 
by X 2 < Xmin + V^data I see text an d Figure[2]for details. (1) S/N ratios near 
the inner worki ng angle are too lo w for optimization of rrj ; we adopt the value 
of 40 AU from Chen et al. 1201 1). (2) The eccentric disk is approximated as 
a symmetric disk offset from the star along the line of nodes. Positive values 
of u indicate a disk center to the West of the star. 



flux to stellar flux of 3.3 [2.9, 3.8] x 10~ 3 . The minimum % 2 
achieved is 1 142 = 1 . lA^ata- The well-fitting family is charac- 
terized graphically in Fig. [3] and summarized in Table [TJ The 
parameters i and g are somewhat degenerate; higher inclina- 
tions mimic enhanced forward scattering. 

In a final step, we explore small non-zero eccentricities of 
the disk while keeping the best-fit values of i, S out , and g 
fixed. We approximate eccentric model disks by translating 
a circularly symmetric disk by an offset u along the line of 
nodes (roughly corresponding to the projected disk's major 
axis), for an eccentricity of e — uja. Positive values of u 
are assigned to offsets to the West. We find a best-fit off- 
set of u — +2.4AU = 20 mas, with a well-fitting interval of 
[— 1.5, +6.3] AU. While the best-fit offset corresponds to an 
eccentricity of e = 0.06, the well-fitting range is consistent 
with zero eccentricity. 

Figure |2] illustrates this best-fit model and its subtraction 
residual. The model fully explains the observed image mor- 
phology down to the noise threshold. 

4.2. Surface brightness profiles 

As a complementary approach, we directly measure the sur- 
face brightness profile of the debris disk. To this end, we con- 
volve the PCA-ADI image with a 5-pixel diameter circular 
aperture and average azimuthally in radial sectors of 9° cen- 
tered on the position angles of the traces of the imaged disk. 
The results are shown in Figure |4&. The profiles are well- 
described by a power law with a slope of 5d a ta = — 2.7 ±0. 1 . 

While the best-fit model disk found in Section 14.11 has an 
intrinsic surface brightness slope of S ou t = — 3.2, the forward- 
modeling of the effects of PCA-ADI data reduction on that 
disk yields a softened slope of S — —2.6 ± 0.3, consistent with 
our data. 
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FIG. 4. — Intensity and polarization profiles of the HIP 79977 debris disk, 
(a) Surface brightness profiles measured in wedges of 9° along the eastern 
and western traces of the debris disk image between position angles [113°, 
122°] and [286°, 295°], respectively. The image was smoothed with a 5-pixel 
diameter circular aperture prior to extraction. The solid curves with dotted er- 
ror sheaths represent the profiles measured on the May 2012 data processed 
with PCA-ADI. The errors were obtained from evaluating a number of 9° 
wedges placed in the background at disk-free position angles. Power-law 
fits to those profiles are overplotted as dashed curves. The fitted slopes are 
—2.6 ±0.1 on both sides of the disk. Furthermore, the surface brightness 
profiles of the best-fit model disk before (top solid curves) and after PCA- 
ADI treatment (bottom dash-dotted curves) are shown for comparison. The 
latter are offset downwards by +1.25 mag/arcsec- to avoid blending with the 
profiles measured from the data. The pure model profiles are not artificially 
offset; the difference of 1-2.5 mag/arcsec 2 with respect to the measured disk 
profile represents the flux loss incurred in PCA-ADI. (b) Degree of scatter- 
ing polarization p at a function of separation, calculated from the polarized 
intensity image (P) measured in the July 2012 data and the full intensity / im- 
age of the best- fit disk model derived from the May 2012 data. The profiles 
are measured in 6° wedges and smoothed with a median filter of 1 1 pixels 
(fs 2 FWHM). The polarized disk flux is lost in the background noise beyond 
l!'5 (cf. FigureGJ). 



4.3. Linear polarization 

We measure the polarized intensity profile in sectors of 6° 
centered on the disk traces in the fully reduced P image, from 
which the RMS halo from residual speckle noise has been 
subtracted. We confirm that the linear polarization is oriented 
tangentially with respect to the star (Fig. [2), as is expected for 
scattering polarization. To calculate the degree of scattering 
polarization, we divide the polarized intensity profile by the 
full intensity profile of the best-fit model disk evaluated in the 
same sectors. 

The results are presented in Figure |4j). We find that the scat- 
tering polarization increases from ~10% at 0"5 to ^45% at 



1"5. This might be due to scattering angles at large separa- 
tions for the considered position angles being constrained to 
^90°, where scattering produces the highest amount of po- 
larization. Close to the star, anisotropic forward-scattering at 
angles > 90° yields a surplus of unpolarized light, diluting the 
degree of polarization. 

4.4. Point sources 

Two faint point sources are clearly visible within 3" of HIP 
79977 in the ADI data, one of which is shown in Fig. [TJl. 
A preliminary astrometric analysis based on marginal detec- 
tions in archival Gemini NICI data identifies both as likely 
unrelated background stars, though more precise astrometry 
is needed for a decisive result. Furthermore, subtraction of 
the best-fit model disk leaves behind a point-like signal of 
4.6(7 significance at a separation of 0"5 (cf. Fig.|2}l). If con- 
firmed, such a signal may represent scattered light from a lo- 
calized clump in the debris disk, or the thermal e mission from 
a you ng planet of 3-5 Jupiter masses (based on Baraf fe et al.l 
2003 and partial self-subtraction correction). While we detect 
no warps or gaps in the disk, this does not preclude the pres- 
ence of a planet of the proposed mas s and separation rang e 



tion range 
aDl2012h . 



orbiting outside the disk's source ring (Thebaultet 

5. DISCUSSION 

The results of our imaging characterization of the de- 
bris disk around HIP 79977 fit very well with the emerg- 
ing standard picture of debris disks, where dust is produced 
in a ring of colliding planetesimals and then distributed in- 
ward and outward by radiation forces. The slope of the 
scattered-light surface brightness is measured to be within 
[—3.6, —2.9], which is consistent with the theoretical value of 
—3.5 typical for a radiation pressure driv en eccentricity dis- 
tribution of particles near blow-out size (iRrivoy et al.l 120061 : 
IStrubbe & Chiangl l2006t iThebault & Augereaul 12001 . No 
significant eccentricity of the disk as a whole is measured, po- 
tentially setting HI P 79977 apart from other disks like t hose 
around Fomalhaut (Chiang et al. 20 09t lAcke et al.ll2.012h and 
HR 4796 A ( Thalma nn et alJl201 lh . However, future observa- 
tions of HIP 79977 resolving the source ring may reveal asym- 
metries below the current uncertainty threshold of e < 0.16. 

The derived Henyey-Greenstein parameter of g w 0.45 
is consistent with m easurements of cometary dust 
dKolokolova et alJ 12004"). Likewise, a maximum degree 
of polarization of 45 ± 15% in //-band agrees well with an 
extrapolation to higher phase angles of obser vations of high- 
polar ization comets like comet Hale-Bopp (Hasegawa et al. 
Il997h . 

As a caveat, we note that other well-fitting disk architec- 
tures may exist outside the limitations of our model assump- 
tions, and that the possible degeneracies of PCA-ADI images 
of disks have not been investigated in depth. 
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